Study of impurities in spin-Peierls systems including lattice relaxation 
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The effects of magnetic and non-magnetic impurities in spin-Peierls systems are investigated allowing 
for lattice relaxation and quantum fluctuations. We show that, in isolated chains, strong bonds form 
next to impurities, leading to the appearance of magneto-elastic solitons. Generically, these solitonic 
excitations do not bind to impurities. However, interchain elastic coupling produces an attractive 
potential at the impurity site which can lead to the formation of bound states. In addition, we 
predict that small enough chain segments do not carry magnetic moments at the ends. 

PACS: 75.10 Jm, 75.40.Mg, 75.50.Ee, 64.70.Kb 



Quasi-one dimensional (ID) spin-Peierls systems at- 
tract intense experimental and theoretical activity for 
their fascinating magnetic properties. Such systems usu- 
ally consist in weakly coupled spin-1/2 Heisenberg chains. 
Due to spin-phonon coupling, these materials undergo at 
low temperature a transition towards a phase exhibit- 
ing a lattice dimerization and a spin gap Q. Inorganic 
compounds like CuGeC>3 are easily doped by magnetic or 
non-magnetic impurities by substituting a fraction of the 
spin-1/2 Cu 2+ ions by spin-0 Zn 2+ or spin-1 Ni 2+ ions. 
As shown by magnetic susceptibility measurements Q 
and inelastic neutron scattering experiments dop- 
ing with impurities leads to a rapid collapse of the spin 
gap. Competition between the spin-Peierls phase and 
a new antiferromagnetic (AF) phase induced by doping 
has been established by magnetic susceptibility measure- 
ments pjfl] , specific heat measurements , neutron scat- 
tering H and NMR experiments Q|. These experiments 
suggest that magnetic moments and enhanced staggered 
spin correlations are induced by impurity doping. 

From a theoretical point of view, the relevant phonons 
in dimerized quasi-lD compounds like CuGe03 or 
NaV2C>5 are often considered as three-dimensional, an 
assumption which, a priori, would justify a classical 
treatment of the lattice. In the dimerized AF Heisen- 
berg chain, a model widely used in the literature to de- 
scribe these materials, one introduces a fixed dimeriza- 
tion J(l ± 5) of the magnetic exchange integral leading 
to the opening of a spin gap oc <5 2 / 3 . In CuGeC>3, it 
is believed that magnetic frustration (i.e. AF coupling 
between next nearest neighbor sites) plays a role. [|]f]] 
The lowest energy excitations of the dimerized Heisen- 
berg chain consist of spinon-spinon bound states [ fb||li"| 
lying below the two-magnon continuum. 

Extensive work on the effect of impurities in dimerized 
spin chains have been carried out Jl2|-[l4j]. The intro- 
duction of vacancies creates finite chains. So far, those 
studies ignore the lattice dynamics which, physically, is 



justified only when the elastic coupling to the neighboring 
chains is large enough. In this finite chain can end 

by either a "weak" or a "strong" bond depending on the 
sign of the dimerization on this bond. These two types 
of boundaries show very different magnetic properties: in 
contrast to the strong bond edge, a weak bond edge can 
localize a S — 1/2 magnetic excitation. This effect is 
responsible for the presence of strong AF correlations in 
the vicinity of weak edges. |ll| 

In the approach discussed above, the effects of impu- 
rities have been considered under the assumption of a 
space- and time-independent lattice dimerization. How- 
ever, due to the magneto-elastic coupling, the presence of 
a spin-1/2 excitation is expected to, locally, distort the 
underlying lattice creating an elastic soliton pq . Such 
effects were recently investigated in the context of the in- 
commensurate phase of spin-Peierls systems under mag- 
netic field |jl|[l8|]. In this Letter, we investigate lattice 
relaxation effects in Heisenberg chains in the vicinity of 
spin-0 or spin-1 impurities by exact diagonalization (ED) 
and quantum Monte Carlo (QMC) simulations. The lat- 
tice is treated either classically in the adiabatic approx- 
imation ||l7f , i.e. allowing for non-uniform static lattice 
distortions, or in a fully quantum mechanical way [jl9 20 
introducing, in addition, the lattice dynamics. In an iso- 
lated chain, we have found that strong bonds form next 
to an impurity. For nonmagnetic impurities, solitonic ex- 
citations do not bind to the impurity. However, in this 
case, the interchain elastic coupling generates an effective 
impurity-soliton attractive potential which leads to the 
formation of a bound state. The spatial extension of this 
bound state is governed by the strength of the interchain 
coupling. 

The model we first consider is purely ID and includes 
a classical lattice distortion. 
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where the second part corresponds to the elastic energy 
lost within the chain. The role of the interchain elastic 
coupling H± will be discussed later. The classical bond 
modulations Si have to be determined from a minimiza- 
tion of the total energy. The pure system is dimerized 
with Si = S(-l)\ 
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FIG. 1. QMC calculation (at T = 0.05) of the modulation 
Si (a) and the local susceptibility Xi (b) (decomposed into its 
uniform Xi an d alternating Xi components) on a L = 79 open 
chain. Data for K = 1 and K = 2.5 are shown. 

We investigate the effect of spin-0 impurities by con- 
sidering open chains. Such chains with an odd number of 
sites must contain at least a spin-1/2 excitation. Previous 
QMC simulations supplemented by a self-consistent de- 
termination of the equilibrium distortion pattern Si (see 
Rcf. jL7) for details) have been extended to this new phys- 
ical situation. Results for a L = 79 sites open chain 
shown in Fig. [l] for various parameters reveal the exis- 
tence of a single solitonic excitation located away from 
the chain edges. The zero temperature local susceptibil- 
ity Xi — J2j {SfSj) at srte * corresponds physically to 
the average value of Sf with respect to the global orienta- 
tion of the total S z spin component. |^] Fig. 0(b) show 
that Xi oscillates rapidly between positive and negative 
values (large staggered component) and has the largest 
amplitude of both its uniform and staggered components 
in the region where the dimer order parameter is sup- 
pressed, i.e. around the soliton. This result is quite 
different to that seen in fixed dimcrization calculations 
p3[ where it was observed that spin-1/2 excitations are 
bound to the chain edge. For increasing K (i.e. for de- 
creasing spin-lattice coupling), the width of the soliton 
increases and the solitonic pattern continuously evolves 
into a sinusoidal distortion as expected in the weak cou- 
pling limit. It should be stressed that different QMC runs 
lead to random degenerate equilibrium solitonic patterns 



centered on different sites in a wide area around chain 
center. However, no changes occur at the two edges of 
the chains which, systematically, end with a strong bond 
(i.e. Si > 0). This clearly indicates that there exists a 
short range repulsion between the impurity (edges) and 
the soliton. One should notice that, independently of the 
presence of the soliton (i.e. for odd or even chains), |<5j| 
increases in the close vicinity of the edges, in contrast 
with the analytical results of Ref . . 
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FIG. 2. Monte Carlo calculation (at T = 0.05) of the 
modulation 5i (•) and the local susceptibility Xi (°) on a 
L = 80 closed chain, K — 2, and a spin-1 impurity at site 
i = 40. 

The case of a spin-1 impurity has also been considered 
by assuming, for simplicity, the same exchange integral 
and elastic constant on the two bonds on each side of the 
impurity. Previous calculations [p2|] assuming a uniform 
dimerization of the chain have shown that spin-1 impu- 
rities lead to more localized states than spin-0 impurities 
(static vacancies). The QMC results shown in Fig. |2| re- 
veal that the two bonds next to the spin-1 impurity be- 
come especially strong indicating that the impurity and 
the two neighboring S=l/2 spins form an effective spin-0 
defect leading qualitatively to the same physics as in the 
case of a S — vacancy. Indeed, the solitonic pattern 
shown in Fig. |^(b) resembles the ones obtained previ- 
ously. However, it should be noticed that the new profile 
is not completely symmetric and that the soliton is al- 
ways located close to the impurity. This signals a small 
attraction in the vicinity of a spin-1 impurity. This might 
be due to the fact that the three-site system formed by 
the S = 1 impurity and its two spin S = 1/2 neighbors 
spends most of the time (i=a 80%) but not all the time in 
the 5 = state. 

In order to investigate the role of the lattice dynamics, 
next we generalize the previous approach by assuming a 
coupling to dynamical phonons, 

H \\ = J^(l+5(^ + ^ t ))S l -S l+1 +fi^ t 6 J . (2) 

i 

For sake of simplicity, optical dispersionless modes are 
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considered here. While the adiabatic treatment discussed 
above is justified in the CI — > limit, the lattice dynam- 
ics can not be neglected when fi and J becomes com- 
parable. In this case, the lattice modulation can be de- 
fined by Si = g(bi + b\\. The treatment of the phononic 
degrees of freedom will relie here on a variational ap- 
proach pj| which gives accurate results [ jO^ . Previous 
calculations |l9| have shown that the lowest S —1/2 ex- 
citations of this model correspond also to massive solitons 
and antisolitons. Furthermore, soliton and antisoliton do 
not bind in the strictly ID case. 
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FIG. 3. Modulation Si (a) and spin density (Sf) (b) in 
the GS of a L — 13 sites open Heisenberg chain coupled to 
dynamical phonons. Parameters are shown on the plot, (c) 
Finite size scaling of the soliton-impurity binding energy in 
the case of a vacancy (as indicated on the plot). 

The results obtained by ED for an open chain with 
an odd number of sites are shown in Fig. |^. The lattice 
modulation pattern and the spatial variation of the spin 
polarization are very similar to the results obtained in the 
adiabatic treatment of the lattice. In particular, strong 
bonds (Si > 0) also form at the chain ends and a soliton 
appears in the middle of the chain. 

In order to get information on the effective interac- 
tion between a soliton and the chain ends, it is instruc- 
tive to define the soliton-impurity binding energy on a 
L = 2p + 1 chain as E B (L) = E IS (L) - E^{L) - e S - ej, 
where Ejs{L) is the ground state energy of an L-site 
chain with an impurity, Eq (L) is the energy of the pure 
system obtained for even number of sites and interpolated 
to L, and es (ej) is the extrapolated soliton (impurity) 
energy (see Ref. ^M). The finite size scaling of Eg is 
shown in Fig. ^| and reveals no binding in the thermody- 
namic limit. The conclusion is also similar when a finite 
magnetic frustration a = J2 /J ^ is considered (in that 
term Tip = J 2 £^ Si ■ S^+2 is added to Hamilto- 
nian (EJ)). In contrast, analogous calculations for S = 1 



impurities indicate a nonzero binding between the impu- 
rity and the soliton. These results including quantum 
lattice fluctuations are consistent with the previous ones 
in the adiabatic approximation. 
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FIG. 4. (a) Local dynamical spin structure factor Su(u) 
calculated on the various sites (indicated on the plot) of a 
L = 11 spin chain coupled to dynamical phonons. (b) Dy- 
namical spin structure factor S(q, to) for the same system as 
(a) (q in units of 2ir/L). Parameters are as in Fig. H. 

The local dynamical spin-spin correlation function is 
shown in Fig. ^(a). The reminiscence of the spin gap of 
the pure chain is clearly seen at an energy u> ~ J. How- 
ever, spectral weight appears at much lower energy. It 
can be attributed to the soliton excitation which behaves 
as a S = 1/2 object weakly connected to the rest of the 
system, consistently with the behavior discussed above 
(Fig. |l|) . If one labels the sites from 1 to L starting from 
the left end of the chain, we observe a large low energy 
weight at the "odd" positions, i = 2k + 1 due to the 
dimerization pattern, with the largest peak at the clos- 
est site to the center. In fact, this can be qualitatively 
understood by assuming that the (free) soliton can move 
by hopping "over" a strong bond. This feature also man- 
ifests itself in the low energy peak in i5(q, lu) near q = tt 
as it can be seen in Fig. ^(b). This is similar to what 
has been observed in fixed dimerization calculations. |l^ | 
The remnants of high energy branch of the pure system 
can be still seen. 

Lastly, we investigate the role of a realistic interchain 
coupling H±. To illustrate the role of H±, let us con- 
sider the physical situation of a finite chain cut by two 
spin-0 impurities at its ends and immersed in the bulk. 
In the dimerized phase, the neighboring chains produce 
a q = 7T potential of the form, H± = K L J2i SiSf xt . Here 
the modulation of the neighboring chains is treated in 
the mean-field approximation, i.e. Sf xt = (— l) l <5o, but 
the full spatial dependence in the chain with the impuri- 
ties is retained. The amplitude of the external potential 
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FIG. 5. (a), (b): QMC calculation (at T = 0.05) of the 
modulation Si vs i on a L — 79 chain with OBC coupled to a 
small external dimerization. (c), (d): the same for a L = 80 
chain. In this case Sf vs i is also shown. 



is then related to the elastic constant K± between the 
chains by K' = K±Sq. In the case of finite chains with 
an odd number of sites, the external potential tends to 
form a weak bond on, let say, the left end. Therefore, the 
soliton will experience a confining force proportional to 
its separation from the left end. This attractive potential 
originates physically from the misfit between the dimer- 
ization pattern on the left side of the soliton and the 
dimerization pattern of the bulk. Our numerical calcula- 
tions shown in Fig. || confirm this intuitive picture. The 
equilibrium position is obtained when the small confin- 
ing potential is equilibrated by the short range repulsive 
potential created by the impurity. 

The case of an open chain with an even number of 
sites is also particularly interesting (Fig. |||(c),(d)). If the 
external potential is out-of-phase with the open chain 
dimerization, the external potential will lead to the for- 
mation of a soliton-antisoliton (ss) pair in the center of 
the chain. For increasing K', the two S = ±1/2 exci- 
tations migrate towards the chain ends forming two lo- 
calized excitations. These calculations again support the 
fact that soliton-impurity bound states are stabilized by 
the interchain elastic coupling. j2J| Because of the finite- 
ness of the energy cost associated with the formation of 
the ss pair, one can, on general grounds, deduce the exis- 
tence of a critical value K' c of K' in such a phenomenon. 
A comparison between the energy cost oc 2es and the 
transversal elastic energy gain cx K'(L — 2F), where T 
is half the soliton width, leads to K' c ~ 2es/(L — 2T). 
Our numerical calculation confirms this prediction. Al- 
ternatively, this implies that, when two impurities are 



sufficiently close to each other, the soliton and the anti- 
soliton bound to each of them can annihilate each other 
leading to the disappearance of the magnetic moments. 
This situation would be more likely to appear for larger 
impurity concentration, i.e. for short chain length on av- 
erage. Note that such features are almost entirely missed 
in fixed dimerization calculations. 

IDRIS (Orsay, France) is acknowledged for allocation 
of CPU-time on the CRAY supercomputers. 

Note added - After completion of this paper we learnt 
of a related calculation on spinon confinement |2q] . 
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